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ABSTRACT 


Contraction theory is applied to an iterative formulation 
of electromagnetic scattering from periodic structures and a 
computational method for insuring convergence is developed. A 
short history of spectral (or k-space) formulation is presented 
with an emphasis on application to periodic surfaces. The 
mathematical background for formulating an iterative equation 
is covered using strai ghtf orward single variable examples 
including an extension to vector spaces. To insure a 
convergent solution of the iterative equation, a process called 
the contraction corrector method is developed. Convergence 
properties of previously presented iterative solutions to 
one-dimensional problems are examined utilizing contraction 
theory and the general conditions for achieving a convergent 
solution are explored. The contraction corrector method is 
then applied to several scattering problems including an 
infinite grating of thin wires with the solution data compared 
to previous works. Problems associated with extending the 
contraction corrector method to two-dimensional iterative 
formulations are outlined including the benefits of applying 
this process to difficult practical problems such as knitted 


mesh surfaces. 
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I. INTRODUCTION 


The electromagnetic scattering properties of periodic 
structures have been of interest to engineers and scien- 
tists for many years and various methods for solving the 
scattering problem have been presented [l] - [8] . These 

methods include solutions generated from static approxima- 
tions, averaged boundry conditions, the method of moments, 
physical optics, the geometrical theory of diffraction, and 
a combination of these techniques. However, certain fre- 
quency regions or certain geometrical configurations (or 
both) cause grave difficulties which cannot be overcome by 
the methods mentioned above .- 

Tsao and Mittra P] presented a novel iterative 
technique based on earlier spectral approaches which 
can be traced to original work performed by Bojarski [10] . 

This technique, called the spec tral - iteration approach, ex- 
tended the scattering solution capability to regions pre- 
viously untouched by other methods. In particular, Tsao 
and Mittra [9] examined the scattering from periodic struc- 
tures and their applications to frequency selective sur- 
faces. The present author's area of interest is the pro- 
blem of determining the reflection propeties of mesh sur- 
faces [ill . Mesh surfaces are used in many applications, 
but the most current application is for conducting reflec- 
tors on space-born antennas. The mesh surface is a complex 



periodic structure whose reflecting properties are not 
readily analyzed by the methods mentioned above. For in- 
stance, an attempt to solve the mesh problem by the method 
of moments would require a special set of basis functions 
as well as an enormous amount of computer memory. Analysis 
of the mesh surface by the spectral -iter ation technique 
would be of great value since the method is essentially in- 
dependent of geometry, i.e., does not require explicit 
knowledge of appropriate basis functions and does not re- 
quire extreme amounts of computer storage. However, the 
basic iterative scheme suffers from convergence problems 
that are associated with most iterative formulations [12] . 
This work treats basic iterative techniques, presents a 
background on the convergence problems associated with 
iterative techniques, details the formulation of a correc- 
tive scheme to iiisure convergence of the iterative tech- 
nique, and applies a version of the corrective scheme to 
the specialized problem of a parallel wire grating. 

Current techniques for solving complex scattering pro- 
blems are limited by the difficulties outlined above. The 
technique described herein provides a basis for the further 
development of solutions to complex problems. Additionally 
the solution methodology provides a useful concept which 
may be applied to problems outside the realm of electromag- 
netic scattering. 



Portions of the background material presented below 
follow directly from the referenced works and credits are 
listed for these. The notation used herein is self con- 
tained with attempts to follow the references as closely as 
possible to maintain a common base, but allowing for differ 
ences to insure clarity of the equations. This is included 
for consistancy, completeness, and continuity. The reader 
is advised to consult the references for greater detail. 
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II. REVIEW OF THE FUNDAMENTAL FORMULATION 

Bojarski [10] considered that field quantity yCx) and 
a source density w(x) were governed by the differential 
equation 

Lj(x)— Lw(x) (1) 

Y w 

where the form of the operators for an n-dimensional 

; p-1,2, — -n (2) 

; p=l , 2 , n (3) 

subject to the constitutive equation 


problem are 


-V /£_ \ 

' T= i?o MVi 


and 




w(x) = q (x) •¥ (x) 


(4) 


The generalized integral representation of equation (1) 
is given by 



x’ ) 


w(x ' )d n x ' + v- __ (x) 
me 


(5) 


where g(x, x' ) is the appropriate Green function satisfying 
equation (1) with w(x) replaced by <$(x) , the Dirac delta 
function and being the externally applied field. 
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Equation (5) is the basic equation which is to be solved by 
application of a transform technique. 

The Fourier transform of a function f(x) is given as 


f (k) = 



f (x)e 


jk*x 


dx 


( 6 ) 


and the transform pairs formed are noted as f(k) «*■*’ f(x). 

Thus, by taking the Fourier transform of equations (1), (2), 
and (3), the k-space formulation of the problem becomes 


L, „(k)?(k)-L (k)w(k) (7) 

y w 

~ m . 

L (k)=£ a.(jk) 1 (8) 

i=o 


L w (k)=£ b (jk) 1 
w i=o 1 


(9) 


where j is the usual imaginary counter-clockwise rotation 
operator. The integral equation (5) then becomes 


y (k)=g(k)w(k) + \ nc (k) 


(10) 


still subject to the conditions of equation (4). This 
approach yields two algebraic equations which may usually 
be solved with much less effort than equation (5). With 



this thought in mind, we now examine the application of 
this procedure to the problem of electromagnetic scattering 
from a periodic structure. 

The electric field E generated from an equivalent mag- 
netic source R can be represented by 


E(x,y,z)= — ^ yxF(x,y,z) 


( 11 ) 


where F is the associated electric vector potential of the 
source and e is the permittivity of the medium [13] . The 
relationship between F and K is established with the use of 
the position vector ? and the free space Green function 


G = 


,-j (k-r) = 

■ I 

4tt Irl 


( 12 ) 


by 


■S' 


F(r ) = / G(r,r ' )*K(r' )dr' 


(13) 


A ^ 

where k and r are the respective unit vectors of the problem. 
From this, the magnetic field intensity H can be derived 
from Maxwell's equations and is given by 


H(x,y ,z)=-joje F(x,y,z) + 


yy.F(x,y ,z) 
jcjy 


(14) 


where p is the permeability of the medium. For the xy 
planar set of magnetic currents shown in Figure 1, F =0 is 
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obtained from equation (13) since z=0 and G is then a 

function of x and y only. Allowing the medium to be that 

of free space, i.e. e= e and y = y , and since F =o, equa- 

oo z 

tion (14) is expanded in Cartesian coordinates x and y to 
yield for z=0 


H s (x.y)= 


3 coy. 


[( 


. a,_ 

k o F * + 


X 

— y, _. 

3 X3 y 


x 

3F 


7 ) 



a- a, 

3F x 8F y \ 

3x3 ^ + 17 - ) 



(15) 


where k Q = is the propagation constant. The quanti- 

ty given by equation (15) is a general representation for 
the scattered magnetic field intensity generated from a 
planar source of magnetic current in the xy plane. 


Consider the planar periodic surface shown in Figure 1 
to be the source distribution for the magnetic field of 
equation (15). Upon substituting equation (13) into equa- 
tion (15) and taking the Fourier transform of equation (15) 
we obtain 


H (a »S )- 

5 


3<uy o mn 


, 2 2 
k -a „ 
o mn 


*’ ct mn^mn 


a mn P 'mn 


? 2 

k -g - 
o mn 


(16) 


^hnn’ B mn j ^ ^mn’ ^mn j 


J(a x 
- mn 


3 y ) 
mn^ 1 



ORIGINAL PAGE - IS 
0(F POOR QUALITY 
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Figure 1. a.) The planar periodic surface indicating 

electric current density J over the 

mn 

aperture portion, 
h.) Unit cell definition. 
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for the transformed magnetic field where 

<W — T ' k o sin0 oos ® 


(17) 


g = cot^-k sinSsinS 

p mn c a o 


(18) 


are the Floquet modes [14] and 


5(a ,3 ) = 

ran mn 




2(k 2 -a 2 -g 2 } 

o mn mn 


(19) 


is the Fourier transform of the dyadic Green function [9] . 
The discrete nature of equation (16) results from taking 
the Fourier transform of a periodic function [15] . The 
final step in completing the process of obtaining a useable 
form of the magnetic field intensity in the aperture is 
accomplished by taking the inverse Fourier transform of 
equation (16), noting K=E Xz, using the equivalence theorem 

a 

and applying the appropriate boundary conditions on 
H S (x,y) at z=0 [9] . This leads to 


H 


tine 


(x,y) 



2 2 
'o ~ a mn 

a mn^'mn 


( 20 ) 


es 

G ( 


a mn 


®mn^ 




8 mn > ( a mn x+ } 
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where tf„ + (a ft % is the transformed electric field in 

a mn, mn; 

the aperture and ^ t ^ nc ^ x »y^ is the incident tangential mag- 
netic field. The inverse of equation (20) is obtained by 
extending the operations over the entire unit cell. The ex 
tended operator is necessary to insure the operations in- 
dicated by equation (20) are in the domain of the solution. 
That is, the boundary conditions of the problem are satis- 
fied. This extension is accomplished by including the elec 
trie current density over the conducting portion of the 
unit cell. The electric current density is included by in- 
troducing the truncation operator and the complement trun- 
cation operator defined .as 

T(f(r)) = 0 for r on the conductor 

( 21 ) 

T(f(r)) = f(r) for ? in the aperture 

and 

T (f(r)) = 0 for r in the aperture 

. _ (22 > 
T (f(?)) = f(r) for r on the conductor 
c 

This allows equation (20) to be extended over the unit 
cell as 



with J(x,y) reprsenting the electric current density on the 
conductor. Equation (23) is a form which can readily lend 
itself to solution by iterative techniques. 
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III. MATHEMATICAL FORMULATION OF ITERATIVE EQUATIONS 
CONTRACTIONS AND FIXED-POINT THEORY 

The use of iterative equations to solve mathematical 
problems has been around for many centuries dating back to 
B.C. 600 when this technique was used to solve problems 
such as determining the square root of three. Beginning 
with this early exploration of iterative techniques came 
the plights of correct problem formulation and convergence 
of the solution. Later methods, such as Newton's method of 
solution, alleviated such problems but often forced other 
constraints on the problem. One such constraint is that 
the initial solution estimate needs to be in a region 
reasonably close to the desired solution. Tsao and Mittra 
[9] used this spirit of a priori knowledge in the form of 
a variational correction to remain near the solution point 
in their problem formulation. However, as will be demon- 
strated later, the cases presented in their work did not re- 
quire such a correction and, indeed, they even noted where 
the iteration method failed even using such a correction. 

The idea of using the iterative scheme to solve gener- 
al scattering problems leads one to seek a way of determin- 
ing the criteria for convergence and a way of generating a 
convergent iterative formulation. Fixed-point theory and 
contractor theory [16, 17] are applied to formulate a con- 
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vergent iterative solution to the scattered fields derived 
in Chapter II. 

The general idea of fixed-point theory revolves around 
finding the solution to an equation such as 

f(x) =y=x^-6x+5= (x-5 ) (x-1 ) = 0 (24) 


This is accomplished by forming an iteration function as 


x n+l “ 8<V <25) 

and finding the intersection y=g(x n > and y=x^ which yields 
the fixed-point and solution to equation (24). This is 
accomplished by starting with an initial guess x q , substi- 
tuting into equation (25) to generate x^ , and continuing 
this process until x n+ i =x n +e where e is the allowed error. 
At this time x=x n+ ^ is the approximate solution to equation 
(25). The conditions for a convergent solution are in 
general : 


a. ) on a closed interval I, g(x) maps I into itself 

b. ) g(x) is continuous on I 


and c.) g(x) is differentiable on I and 


dg(x) 


dx 


< 1 


Thus the following theorem may be stated 

Theorem 1 . Let g(x) satisfy conditions a.), b.), and 

c.). Then g(x) has exactly one fixed 
point x* in I and starting with any x Q in 
I, the sequence x Q , x^, ... , x n generated 
by equation (25) converges to x*. 
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The process outlined by the above theorem is illustrated in 
Figure 2. The fixed-point x* is often said to be a point 
of attraction for g(x). x n+ i is also said to contract to 
x*. 


To illustrate the use of this theory, one root of equa- 
tion (24) will be determined. Let 


x 2 +5 

x n + rS (x n>= -V- 


(26) 


and find 


d 

Hx 


g(x) 



(27) 


then, the interval I is defined as x«I and IcC-3,3] . Note 


that the known root x=l is in this interval and 
d 


=1/3. We then expect that the iteration function 
x=l 


35 S (x) 

will converge to the fixed-point x*=l which is the correct 
solution. The iterative process is contained in Table 1. 
The question then arises, what happens when the sequence 
will not converge? 


To generate an iterative equation that will convege we 
must meet the conditions listed in Theorem 1 above. As an 
example, let us examine the roots of 

f(x) = - x - 6 = (x+2) (x-3) 


(28) 



ITERATION EXAMPLE 


15 



C XU 


Figure 2. Example of convergent iterative equation g(x). 
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f(x) = x - 6x +5 = 0 


g(x) = x = 


x 2 + 5 


g(x) 


x 0 

2.0 

X 1 

1.5 

x 2 

1.208 

x 3 

1.077 

x 4 

1.026 

x 5 

1.009 

x 6 

• 

1.003 

• 

• 

x n-l 


x*=x n 

1.000 


Table 1. Example of a convergent iterative equation 
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which can be seen to be -2 and 3. 
equation 

« (x n ) ' x n + r x n 2 ‘ 6 
and noting 

|a*e<*>| HHk-2,3’ 1 

we see that the iterative equation will not converge. How- 
ever, by properly forming a new iterative equation G(x) 

that has ^ G(x) x= _2 < ^ * or Hx x=3" <: '*' ’ t '^ le root:s 

may be obtained. Although manipulation of equation (29) 
could yield a suitable G(x) [12] we seek a G(x) that can 
be formed by a method known as relaxation. The relaxation 
process generally uses a portion of the n t ^ 1 iteration along 
with a portion of the iterated n L (i.e. n+1) iteration to 
"relax" the process. The form of the relaxed iterative 
equation is 

G(x)=Rx+(l-R)g(x) (31) 

n n ° n 

where R is called the relaxation constant. The condition 

=0 (32) 

x=x 

n 

may be used to determine an optimum R such that convergence 


k g(x > 


By forming the iterative 


(29) 


(30) 
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of G(x) is assured. Performing the operation of equation 
(32) on equation (31) yields 


dx 


4- G(x)|= R+(l-R) 4- g< x > 


dx 


X=x_ 


= 0 


(33) 


n 


x=x 


n 


and solving for the linear correction R 


R- K ' (x) 

R - rnm 


(34) 


where g'(x)= ^ g(x)l. This process of finding the optimum 

lx=x 

R is called the contraction corrector (a name to be ex- 
plained later in this section) and is used to iteratively 
solve for the root x=3. This is illustrated in Table 2. 
Equivalence between the contraction corrector and Newton's 
method can be shown and is illustrated in Appendix C. 

Note that g(x^) would not converge by itself and thus G(x^) 
is necessary for a correct solution to be obtained. As a 
further note on the iterative process, the conditions that 
g(x) maps I into itself is not trivial. For instance, note 
that 


f(x) = x^ + 2 = 0 


(35) 


will not yield a solution from an interative process unless 
x is allowed to be complex and g(x) or G(x) maps into the 
complex plane. The iterative equation G(x) may even re- 
quire a complex relaxation or complex contraction correc- 
tion constant and precautions governing complex algebra 
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(K n 2 -6) 


n 

X 

n 

R 

n 

g(X n ) 

G (x n ) 

0 

1.0 

2.0 

-5.0 

i 

7.0 

1 

7.0 

1.07 

43.0 

4.24 

2 

4.24 

1.13 

11.98 

3.23 

3 

3.23 

1.18 

4.43 

3.01 

4 

3.01 

1.2 

3.06 

3.0 


n 

*h 

g(x n ) 

1 

3.0001 

3.0006 

2 

3.0006 

3.0036 

3 

3.0036 

3.0216 

4 

3.0216 

3.130 

5 

3.130 

3.798 

6 

3.798 

8.424 


Table 2. Example of a non-convergent iterative equation 
with convergence achieved by applying the con- 
traction correction method. 


g< V ' x n+l " x n ' 6 


n+ 


1 = G(x n ) = Rx n +(1-R) g(x n ) = Rx n =(l-R) 

2 x_ 


R = 


n 


2V 1 
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and operations must be observed. In particular, the ana- 
lytic nature of f(x), g(x), and G(x) in the interval I 
should meet the requirements of Theorem 1. An illustration 
of the nature of these problems is given in Table 3. In 
two cases, a real x Q never maps x n into the complex plane 
and, as such, never yields a solution. Even applying the 

contraction corrector method with a real x does not give a 

o 

solution. Only the one case using a complex x q and the con- 
traction corrector method converges. This follows from the 
conditions of Theorem 1. Since problems have been observed 
in the simple problems above, a proper choice of action is 
to pursue the corrective method in a more general sense. 

Stakgold [l8l has formed a very sound collection of 
theorems and definitions on the idea of metric spaces and 
their transformations, and the following material is con- 
densed from his book with the advice that the reader con- 
sult his work for proofs and greater detail. 

Definition 1. A transformation L of a metric space 

x into itself is Lipschitz continuous 
if there exists a p, independent of u 
and v, such that 

d(Lu,Lv)_< pd(u,v) for all u,v£x 
where d(S,£) is a proper metric in x. When the definition 
above holds for some fixed \p| <1 , the operation of L is 
called a contraction on X. 



X + 2 


1.00 

3.00 

11.00 123.00 

- diverges 


33 

-9.33 0.9 

-3.83 
69 


1.11 3.24 I 2.231 1.81 


4.24 + J4.24 1.88 + J35.95 
2.18 + J2.34 1.81 ♦ J10.24 
1.13 + jl.45 1.20 ♦ J3.29 
0.64 + jl .25 0.85 ♦ jl.59 


1.06 - JO. 07 
1.09 - JO. 13 
1.13 - j0.28 
1.05 - JO. 40 


2.18 + J2.34 
1.13 + Jl .45 
0.64 ♦ jl .25 
0.50 > jl .31 


X* - 0.5 J1.32 

Example of a complex valued iterative equation 
and the use of a complex contraction corrector 
















Theorem 2. 
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Let the operator L be a contraction 
on a complete metric space X. Then, 
u=Lu has one and only one solution 
in X which may be obtained from any 
initial starting point u Q in X. 

Thus, u^ - *" u* (the fixed-point) as 
n-*-» and u* exists and is unique in 
X. 


Note that u may have multiple components, i.e., u=u(z^,Z 2 ».. 
. »z n ) where z^ are complex coordinates. P as given by defi- 
nition 1 is the description of the measure of a general der- 
ivative of an operator. Thus, the process used in the sing- 
le variable examples earlier may be extended to higher di- 
mensional spaces with more general operators. This theorem 
allows that the relaxed corrective iterative process out- 
lined earlier is indeed a contraction, hence the name con- 
traction corrector. This nomenclature is used to differen- 
tiate this scheme from other corrective schemes applied to 
iterative processes. The use of this definition and theor- 
em as applicable to the spectral -iteration approach of equa- 
tion (23) is detailed in the following section. 
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IV. THE PERIODIC STRUCTURE: A 

CANONICAL CASE AND THE REGIONS OF SOLVABILITY 

The total electric field for the structure shown in 
Figure 1 may be expressed in terms of the incident and 
scattered fields as [9 ] 


(v“)* 

' ^mc- 


l 

mn 


a 6 
mn mn 


k 2 a 4 
o - mn 


-k +6 -ai 
o mn nin mn 


G ( a mn» ^mn^ 


( 5 »„„) 'C--* 

a y. 


®mn y ) k z 2 for * >0 

' mn 


mn 




l 

mn 


a 6 
mn mn 

-k 2 + 6 2 

o mn 


k 2 - a 4 
o mn 

-a 6 
mn mn 


(36) 


ra . B \ 

1 * mrj I 


G( mn* mn 



for z<0 


where 


mn 



B *)* 

for k c 

+ Bmn* * 

t K 
0 

N 

*** 

for k Q 


2 > a 2 + B. 2 


mn 


mn 


( 37 ) 


mn 


2 + 6 2 
mn 



/ 
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and E and E represent the transformed electric field 
and are respectively the reflection and transmission coef- 
ficients of the Floquet modes. Deriving H from equation 
(36) and by enforcing the boundary conditions on the tan- 

i 

gential magnetic fields across the cell at z=0 , equation 
(23) is generated again. This exercise is presented to 
complete the field description and to illustrate the rela- 
tionship between the aperture fields and the reflection and 
transmission coefficients as the coefficients will be of 
interest later in the section. Equation (23) is the begin- 
ning point to cast the problem in an iterative form. 

The summation in equation (23) represents a discrete 
Fourier series (DFS) for an infinite duration (i.e., peri- 
odic) sequence [19] . This representation allows for direct 

transformation between the (x,y,z) domain and the (k ,k >k ) 

x y z 

domain. Since the functions are represented by sequences 
with complex exponentials having a periodicity of 27r/m and 
2n/n, one period of the aperture distributions, that is one 
cell of the structure, can be used to completely specify the 
transform. The use of one period to represent a periodic 
function in this manner is known as the discrete Fourier 
transform (DFT). It is important to note that the trans- 
formation mentioned above is exact and this infers that no 
aliasing will occur when the DFT is performed. Aliasing 
can occur when a non-periodic function is truncated and 
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and transformed with a DFT. Since information about the 
original function is lost in the truncation, performing the 
inverse transform can never return the original function. 
However, when using exactly one period of the periodic func- 
tion, the original function can be reconstructed when the 
DFT and inverse DFT are performed. The most common form of 
DFT algorithm is the fast Fourier transform (FFT) in which 
special properties of the DFT are exploited to decrease the 
computation time of the transform. Throughout this section, 
the DFT and its inverse are represented by F and F ^ re- 
spectively. Implicit functional dependence is also used to 
present a "cleaner" form of the equations. 


With the above thoughts in mind, we may solve equation 
(23) in terms of E t as: 






V J) 


H 


tine 


(38) 


and, for clairity write equation (23) as: 


T (J) = T 
c c 


H. . 
tine 



( — 
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G i F 

T(E t ) 


( I 


' . 


(39) 


Note that the subscript indicating aperture field has been 
removed from equations (38) and (39) since the truncation 
operator is included and the tangential field E t over the 
entire cell is used. Solving for E t in equations (38) with 
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T c (J) substituted from equation (39) yields: 



This is the form of iterative equation given by Tsao and 
Hittra [9 ] and is the basic equation to which a scheme 
called the contraction corrector will be applied. The con- 
traction factor must be less than 1.0 to have a convergent 
solution [18] and the contraction corrector scheme will 
allow this condition to be met. 


The iterative equation for one-dimension (E x or E^ 

will converge if the conditions for a contraction given by 
Stakgold [l8l are met. The proper metric for the space X 
and the operators in equation (40) is: 


d(u, v) 




where u = ( K\ > * '* * * anc * v = ^l’ ?2’’" ’ ^ are t * ie 

vectors in X and x, represent the components of each vec- 
tor. The u and v used in the metric d(u,v) may be any u,v 
in X since the contraction must hold for all u,v in X. One 
logical choice for u and v is E^ n ^ and E^ n+ ‘' - \ This choice 
would allow the contraction process to be observed and 
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would give the relative error between each iteration. This 
also allows the regions of convergence to be determined. 

That is, for decreasing error, the solution is converging 
and for increasing error diverging. Since we are interest- 
ed in the contraction to a fixed-point, a wiser choice for 
u and v is and + e) where e is some small complex 

number. Thus, at each component of E^ n ^ we essentially 
have the ability to determine the general derivative of the 
operator at each element. The contraction corrector is a 
vector with one element for each component of E^ n ^' 

The form of is exactly that of equation (31) for a one- 

variable problem extended to an N component one-dimensional 
problem. This powerful concept allows the regions of con- 
vergence to be determined and allows the formulation of a 
convergent solution as 

Tg(n+D = g-(n) . g(n) + ( f _ R< n >) » L £ (n) (42) 

with L being the operator defined by equation (40). Unique- 
ness of the solution is guaranteed since the problem was 
formulated from Maxwell's equations via a vector magnetic 
potential and the appropriate boundary conditions were 
applied. Note that this representation is for a one-dimen- 
sional scattering problem. Additional conditions are nec- 
essary for the general solution of two-dimensional problems 
[12] . These conditions and other considerations for two- 
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dimensional problems are discussed in Section VI. 

The contraction corrector R = 

be determined in a general operator sense to behave as R in 
equation (34). This implies that closer examination of De- 
finition 1 is warranted. The determination of an operator 
derivative must be formulated in terms of a contraction. 

The derivative must be expressed in an analytical form that 
allows the contraction to be determined numerically for 
each R^ of R. Thus, equation (42) will guarantee a conver- 
gent solution when R is properly determined. The contrac- 
tion is related to the .continuity of the operator (trans- 
formation) L through Definition 1 as: 

d(Lu, Lv) ^ P d(u,v) (43) 

with d(£,C ) given by equation (41). The determination of 
the contraction of L at u = u and v =(u + Au) is a represen- 
tation that is thought of in operational calculus as the de- 
rivative of the operator L for vanishing Au. Equation (43) 
is then written as: 

d[Lu, L(u + Au)] £ p d[u, u + Au] (44) 

If the vector u has. m components and equation (44) is expanded 
using equation (41), the result is: 


Ri , R, 


R 


n 


must 
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LtUj+AuJ-LUj 


+ L(u 2 +Au) - Lu 2 


+• • •+ L(u +Au)-Lu 
1 m m 


Jti^+Aui-Uj 2 + j(u 2 +A u )-u 2 2 +*••+! 


(u +Au)-u 

m m 


(45) 


Eliminating terms in the denominator of equation (45) yields 


1 


~ 

Ltu^+AiO-LUj 

2 + L(u 2 +Au)rLu 2 2 ♦•••+ 

L(u +Au)-Lu 2 
m m | 


< p (46) 


Au 


Now, since m>l, a stronger condition is written 

|l(u^+Au)-Lu 1 | 2 + |l(u 2 +Au)-Lu 2 | 2 ♦•••+|L(u m +Au)-Lu m | 

" TauT — < p (47) 

Consider the following: given 
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Is the inequality true? Squaring both sides yields: 
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and the inequality is seen to be valid for non-zero a^. 
Equation (47) is then written as: 


L(u, +Au)-Lu. 


L(u 2 + Au)-Lu 2 


+ • • • + 


L(u +Am)-Lu 
m 


< p (50) 


Au 


and then may obviously be further reduced for each component 
t0 ’ L(u i +Au)-Lu. 

<" A 

(51) 


1 p i 


Au 


m 


noting that .2 P, < p . 

’ i=0 l 

From equation (51) it is seen that a complex-valued p may 

c i 

be defined and the final form of equation (51) becomes: 


L(u.+Au )-Lu. 
i i 


Au 


= p c. 

l 


(52) 


This form is the measure of the derivative of the 
operator L at each component and may be used to solve for 
the contraction corrector R at each component u^ e u. 
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Even though a powerful tool has been developed to 
allow the computation of a convergent one-dimensional solu- 
tion to electromagnetic scattering, the background problems 
surrounding the contraction corrector scheme should be ex- 
amined. This . background material is of interest since the 
breakdown of prior methods and the reasons for non-conver- 
gent solutions gives insight into what particulars of the 
scattering problem cause the difficulties. In particular, 
known geometrical problem areas are: conductor width, cell 

size, aperture size, and the angle of incidence and polari- 
zation of the impinging plane wave. The understanding of 
the difficulties arising for certain values of those para- 
meters can provide vital insight when attempting to extend 
the solution to two-dimensinal configurations, especially 
when mesh structures are involved. The reference point for 
this background investigation is the work of Tsao and Mittra 
[9] . A complete iterative flow is detailed in Appendix B. 

Consider the case of the free-standing strips shown in 
Figure 3. This planar configuration of thin, perfect con- 
ductors was the prime example used in the work of Tsao and 
Mittra [9] . The case of a cell width of 1.4 a. , aperture 
width of 0.6x , and incidence angle 0=0* was of particular 
interest. Although the iteration equation (40), would per- 
form adequately, Tsao and Mittra decided to apply a correc- 
tion step to the iterative process. The idea of being 
"near" the solution in numerical methods (i.e., Newton's 
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Method) to either insure or speed convergence is commom 
place. The correction step was chosen to be a form of the 
Method of Moments [28] . 

The amount of amplitude error of E^ n+ ^ may be calcu- 
lated from processing with equation (40) and using an 

amplitude correction C determined via a variational formula 
calculated with Galerkin's method. This correction is 
given by [9] : 


C 


< 


g< n ) _jji 

tan, tan 


> 


l< W ( n ) 
ran. 


,-l 


E< n) 

tan 


(53) 


where : 
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x^ X X 2 


ds 


aperture 


(54) 


with denoting conjugation and d¥ the outward surface 

normal. The value of C yields a weighted average error 

based on the first moment's variation of the operator using 

E^ n ^ as the test function. C will be equal to one when E^ n ^ 
represents the averaged approximate solution to the problem. 

This statement implies that oscillations of the solution 

could average out and yield a solution that is incorrect. 
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Numerical difficulties could then arise if a problem gener- 
ated such a solution from the iterative process. Great 
care is then warranted with this type of numerical method 
as is true with any numerical scheme. Additionally, as pre- 
viously noted for iterative processes, just being near the 
correct solution in amplitude does not guarantee a conver- 
gent process. 

The geometry given in Figure 3 for a=1.4X, b=0.6A, 0=0° 
and H-wave parallel to strip edges) incidence is a 

problem that converges by equation (40) without the aid of 
the variational correction of equation (53). What then is 
the useful purpose of the correction C? The most useful 
purpose of C is to speed convergence since an iterative 
equation which is a contraction will converge from any ini- 
tial estimate and will diverge if not a contraction. Anoth- 
er useful (but applied with caution) purpose is the indica- 
tion of accuracy of the iterated solution. The problem 
above is taken from Tsao and Mittra [9] . The convergence 
of the problem is shown in Figure 4. Figure 5 illustrates 
the convergence of of the problem using equation (40) 

corrected by C of equation (53). The iterations of both 
Figure 4 and Figure 5 are quite similar; however, the 
shown in Figure 5 are converging more rapidly then the E^ n ^ 
in Figure 4 by one iteration. This is a very small difference 
in convergence rates and indicates that the methods behave 
similarly. 
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Equation (42) is also used to solve the problem above. 
The contraction corrector R is obtained using equation (52) 
along with equation (34). The application of this contrac- 
tion corrector yields results for shown in Figure 6. 

Note that the first iteration is the same as the first 
iteration in Figure 4. However, the second and third itera- 
tions shown in Figure 6 are improved over those in Figure 5. 
Even though this convergence improvement is by one itera- 
tion, it is pointed out that the contraction corrector 
guarantees convergence (except for possible numerical compu- 
tational problems). The case for using the contraction cor- 
rector over the other methods is not well presented in the 
preceeding example. The basic iterative equation and the 
variational corrected iterative equation break down and di- 
verge for various cases. The following material is present- 
ed to establish bounds for the regions of solvability for 
the preceeding problem and to demonstate the usefullness of 
contraction theory in this matter. 

The basic iterative formulation is found lacking for 
various cases. Conductor size, aperture size, cell size 
and incidence angle are the parameters that are varied and 
may cause convergence difficulties. The first case to be 
examined is that of constant Floquet cell width and varying 
aperture size. In paticular, let a = 1.4x and 0.005X £ b 
<_ 1.395X for perfectly conducting strips, incident angle 0=0° 
and H-wave polarization. The contraction factor for this 
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Figure 7. Contraction factor for a cell width of 1.4X, incidence 
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Figure 8, Contraction factor for a strip width of 0.6A, incidence 
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case is graphed in Figure 7. Recall that the solution will 
diverge for a contraction factor greater than one. These 
non-convergent regions are where the contraction corrector 
must be applied to insure covnergence. Next, the cell size 
is varied as 0.9\ <_ a <_ 1.4x for constant conductor size, 
i.e., (a-b) = constant, and 0= 0° with H-wave polarization. 
Figure 8 shows the contraction factor for a conductor size- 
0.8 A . This last example contained a conducting strip that 
was a large portion of a wavelength. The case of a conduct- 
or size 0.1X and cell width 0.2A < a 0.9A for H-wave po- 
larization and 0= 0° is presented in Figure 9. This example 
represents a small strip and the contraction factor associ- 
ated with its geometry. The last case to be examined is co- 
cerned with the angle of incidence of the impringing plane 
wave. A cell width a = 0.9 A , conductor width b = 0.4 A 
and aperture width 0.5 A is chosen for the convergence sur- 
vey illustrated in Figure 10. Note that the contraction 
factor bumps above 1.0 around 0=8° and diverges from 1.0 
for 0>6O°. This indicates that the basic iterative scheme 
will have problems with convergence around 0=8° and will 
not converge for 0>6O°. The above examples offer several 
cases that point out problem areas in constructing a conver- 
gent iterative solution. 

The contraction corrector may be applied to construct 
a convergent solution when the basic iterative scheme 


42 



Q) 

O 

C 

0) 

T3 

■H 

O 

C 

•H 


iH 

o 


o 

II 

e> 


<H 

o 

x: 

■p 

T3 

•H 

5 


ft 

•H 

Vi 

-P 

CO 

OJ 

Vi 

O 

P 

Vi 

O 

-P 

O 

cij 

P 


c 

o 


o 

ClJ 

V 

-p 

c 

o 

u 


O 

a> 

p 

3 

bO 

•H 

ft 


UOJ^30JV»3 


angle of 0=0 with Il-wave polarization and 



Figure id. Convergence survey using contraction criteria 
for 0° <0<90° , 0=O°,a=O.9A, and b=0.4X. 
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breaks down. Consider the case of Figure 3 with a = 1.4 X, 
b = 0.6X and 0= 70° . At this incidence angle, the itera- 
tion equation (40) and variational correction of Tsao and 
Mittra [9] break down and do not converge. Figure 11 illu- 
strates that after six and seven iterations respectively, 
neither method converges. However, Figure 12 reveals that 
the contraction corrector method converges after six itera- 
tions and has a small error over the conductor. Herein 
lies the reason the contraction corrector method is pre- 
ferred. Even when the basic method and variational cor- 
rection schemes fail to converge, the contraction cor- 
rector method assures that solution can be reached and does 
so with reasonable ease. 

The foregoing examples have demonstrated the ability 
of the contraction criteria to be applied to one-dimensional 
scattering problems to determine the regions of solva- 
bility. Additionally, the contraction corrector scheme has 
been utilized to achieve convergence when the basic itera- 
tive schemes failed. The examples presented give indica- 
tions that the contraction corrector can be successfully 
applied to planar wire surfaces. 


I 






ORIGINAL PAGE IS 

OF POOR QUALITY 


46 



PI*!J 3 


Figure lib, Non-convergent results for 
equation with 0=70° . 
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V. THE PROBLEM OF AN INFINITE GRATING 
OF THIN WIRES 


The structure shown in Figure 13 has been studied over the 
years by various techniques to determine the reflection co- 
efficient from the grid [1, 2, 20 ] . The reflection coeffi 
cient may be obtained from the transform of the aperture E 
field in the form of the Floquet modes. For the case of 
cell widths less than A/2, this coefficient is obtained 

from the n = m = 0 element of the E~ variable for the 

v 

- mn 

case of plane wave incidence as [9] 


r 


~ ( n) 
E 9oo 


- (n) 
= E + 
yoo 


E . 
ymc 


(55) 


The first atempt to solve the problem with a A/4 cell and 
A/600 radius wire was performed with the variational correc- 
tion iterative process using E^ n ^ as a test function in the 
variational scheme. This method fails to converge for vari 
ous sampling rates and various numbers of Floquet modes, 
and a more suitable test function is sought [21] . The 
Dirac delta function is chosen as the test function and 
gives acceptable results for small angles of incidence. 

The delta function samples the response and allows point 
matching to be utilized to improve the variational correc- 
tion factor. The convergence properties achieved with these 
two test functions are shown in Figure 14. Since the iter- 
ation method should be convergent for all incidence angles. 
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the development of a different form of correction is a 
logical path to follow. 

The iterative process of equation (40) is one derived 

for E from itself. That is, the E^ n+ ^ was determined from 

E^ n ^ with an intermediate step in the iterative process 

used to determine the surface current density J^ n ^. The 

iterative process then has two distinct parts and it is 

reasoned that two distinct correction factors are required. 

By weighting previous iterations of both E^ n ^ and J^ n ^’ 
is found that convergence can be obtained for small cell 

size and near normal angles of incidence. However, it is 
also determined that the stability and speed of convergence 
is critically dependent on the weights chosen. It is then 
realized that a solution may only be obtained by an "acci- 
dent" as mentioned before. The hit-and-miss approach re- 
quired by the weighted iterations may allow the determina- 
tion of a proper solution but the recursive time required 
could be excessive. These facts then drive the iteration 
equation toward a solution obtained by the contraction cor- 
rector method. 

The cases chosen for examination are x/8, A/4, and A/2 
cell widths with both TE and TM incidence. The TE case cor- 
responds to incidence angles $ = 0° and 0° _< 0_< 90 u while 
the TM case corresponds to incidence angles <5= 90° and 0 U 

< 9< 90° , with either E. or H. respectively co- 

- - inc inc r 7 
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polarized Co the wires. Since the determination of reflec- 
tion from two-dimensional mesh reflector antenna surfaces 
is the application toward which this work is directed, the 
reflection coefficients must be carefully scrutinized to 
ascertain the best definition of reflection. Wait [2] de- 
fines a suitable normal reflection coefficient and a com- 
parison to his results will be made. The specular reflec- 
tion is of interest in antenna work and the construction 
of the geometry required is given in Figure 15 [22] . When 
E. is entirely co-polarized with the grid wires, the two 
reflection coefficients agree in definition. With TM polar 
ization (vertical) incident, a pseudo-Brewster angle is ex- 
pected at certain incidence angles, Variations of magni- 
tude are also expected to occur with changing wire conduct- 
ance. Before examining various cases and comparing the re- 
sults to those of previous authors, the wires must be con- 
verted into equivalent strips to yield a structure equiva- 
lent to that used by the previous authors Cl, 2, 20, 27 ] . 

Harrington [29 ] relates that equivalence between thin 
metallic rectangular strips and small radius wires is ob- 
tained if 

b - 4w (56) 

where b = radius of the wires and w = width of the strips. 
Thus, using equation (56), we are able to use the strip ana 
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lysis to determine the reflection from a grid of cylindrical 
wires. A further characterization of the wires is required 
when the wires are of finite conductivity. Schelkunoff [23] 
gives the internal impedance of a thin cylindrical wire as 


Z 

1 


t 


nm K o ( ^ b) 

2^bK^yb) 


(57) 


where n = in tn/( o + Jcue ) ^ and Y = ip w(c + ioj£ ) % 
r| m °m J m m J m 

with v ,e > and c being the electrical properties of the 
m m ° r r 

wire and K q , being modified Bessel functions. These 
principles may now be applied to the cases detailed above 


The cases of cell widths of X / 8 , A/4, and A/2 are now 
presented with a wire radius of a/ 600. Good agreement is 
obtained between the results of the contraction correction 
method and the results of Wait [ 2 ] for incidence angles of 
both $= 0° and $= 90°. The case of $= 0 U is illustrated in 
Figure 16. Recalling the other definition of reflection co- 
efficient given in Figure 15, the reflection coefficients 
for the case of $ = 90° are shown in Figure 17. The expect- 
ed pseudo-Brews"ter angle is observed as a reflection mini- 
mum near the angle of 67*. The internal impedance of the 
wire becomes important in the problem when the conductivity 
of the wires becomes finite. Instead of the tangential 
electric fields cancelling over the conducting portion of 
the cell, they must combine to support the current at the 
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Figure 16 . Comparison of results obtained from Wait C2} 
and equation (30) for various cell widths 
with P=0° , 0° <0 £90°, and o= 00 S. 
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surface as 

E. . + E. _ = I Z. + 

tine tscat o mt (58) 

where E_ _ = scattered electric field and T = current 
tscat o 

along the conductor. The results of a finite conductivity 
are shown in Figure 18 where equation (58) has been included 
in the iteration equation. The results resemble that given 
by Jordan and Balmain [30] for a lossy ground and follow 
the idea that the grid represents a shunt impedance. Over- 
all, good agreement is obtained with previously computed and 
measured results. 

As a final note, some concern has been expressed over 
the non-convergence of the iterative technique and over the 
fact that increasing the number of Floquet modes not not in- 
crease the accuracy of the solution [24] . The first con- 
cern has been addressed in the preceeding sections while the 
second will be addressed now. When an iterative equation is 
in a non-convergent region, the solution may only be arriv- 
ed upon by an accident of computation unless special methods 
of computation have been included. 

The consequence of adding additional Floquet modes is 
adding additional points that are in the non-convergent re- 
gion. The result of this consequence can range from produ- 
cing a large variation in the iterated solution to helping 




Figure 18. Several cases of finite conductivity with 
0=0°, 0° <©<90° , a=0.25X and b=0.247X. 
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the iterated solution. The control over the results is nil 
with a non-convergent iterative scheme. However, when the 
contraction corrector method is applied, the expected in- 
crease in accuracy with increasing Floquet modes is obtain- 
ed. Of course, a limit on accuracy still exists due to the 
finite register length of the computing machine used. The 
case of a X/4 cell width, X/600 radius wire with infinite 
conductivity and $=0=0° is chosen to illustrate the above 
comments. A plot of converged accuracy is shown in Figure 
19. The usual leveling off accuracy is noted. This is due 
to the higher order Floquet modes having less and less 
effect on the solution, i.e., the farther away wires con- 
tribute less and less to the fundamental cell when a con- 
vergent process is used. 
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Figure 19* Illustration of converged accuracy .for various 
number of terms with a=0.25A, b=0.247A, 0=9=0° 
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VI. EXTENSION TO TWO-DIMENSIONAL PROBLEMS 

The present analysis has been limited to fields having 
only one spatially varying component. When a second spa- 
tially varying component must be added, the difficulty in 
applying the contraction corrector scheme is increased. 

Each field component must have its own iteration equation 
as can be seen from equation (23). These equations are 
coupled, however, and this leads to new constraints that 
must be applied to the problem Cl2l . For iterative equa- 
tions of the form 

x n+l = f(x n'3'n ) 

(59) 

^n+1 * 8 (x n-y n ) 

in a space X, the sufficient but not necessary conditions 
for convergence are 


1 3 f 

+ 

3 f | 

I 3x 


3y | 

1 3g 

+ 

3g I 



3y 1 


P 


< 


P 


62 


for all x,y in X where p< 1. As was observed in the one- 
dimensional case, the conditions for convergence revolve 
around the derivatives of the iterative equations. The de- 
termination of the partial derivative follows the form of 
equation (52) and may be written for f(x n ,y n ) as 


f = 
x 


3 f 
3x 


■V*n'V = 


t(x n + ;ax, y„>-f( yy„> 

Ax 


(61) 


where f indicates the partial derivative of f(x,y) is 
taken. Again, it is seen that contraction theory allows 
the measure of the derivative to be calculated. The condi- 
tions to insure convergence must now be determined. 

The conditions given in equation (60) are sufficient 
and not necessary and as such are fairly strong conditions. 
Thus if a stronger condition is applied, the stronger condi- 
tion is also sufficient. Such conditions are 



These conditions can not be met unless a new set of equa- 
tions are formed using the contraction corrector method. 

The method for the one-dimensional case required the com- 
bining of x n and x r+ j together to form a new iteration equa- 
tion. The two-dimensional case requires the combining of 
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the x and x , along with f , f , g , and g to meet the 
n n T x x y x y 

requirements of equation (60). A form of equation that 
conforms to the plan outlined above for x n+ ^ is 

x n.l' F(x n' y n >=ax n +<1 -“ 6 <63 > 

The first term of equation (63) is that of the one-dimen- 
sional case. The second term of equation (63) is seen to 
be equal to zero when y n approaches the fixed-point y* . 
Similarly, the iterative equation for y n+ ^ is: 

y n + l = G(x n-5 , n ) = tYy n +<l - Y>g<VV ]+[3x n' 3 

(64) 

The complex constants a, 6 ,Y , and 3- are the contraction 
corrector terms that must be utilized to force equations 
(63) and (64) to meet the convergence requirements of equa- 
tion (60). 

The partial derivatives of equation (63) are set equal 
to zero and are found to be : 


> £ X <VV - S 8 X (x n-Y n > * 0 


VW = <1 -“ )£ y (x n’Yn >+B - 6 8y <X n’ y n ) = 0 


From these equations , 8 may be solved in terms of a as: 


(65) 

( 66 ) 
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6 = 


° + (1 - a,f x (x n>V 


VVV 


for g x <x n .y n > + 0 


(67) 


which may be substituted into equation (65) to find: 


a = 


f x ( x *y ) 

CVW - g^(x",y") - WV 

C f x (x n’5'n ) ' ^ ^ 


( 68 ) 


8?V7 


- V vV 


If g (x ,y ) = 0, a is found from equation (65) as 
x n n 


a 


f (x ,y ) 
x v n °n - 


WV - 1 


(69) 


which is seen to be the same form as equation (34) in the 

one-dimensional case. This is expected since F (x ,y ) is 

x n ,J n' 

not affected by g (x , y ). 3 is found from equation (66) 

a n n 

as : 


6 


( a -1) 

1 ~ g 


f 


y 


z 


for ®x (x n , ^n ) = 0 


( 70 ) 


The case g y <x n ,y n ) = 


1 is excluded for the problems of 
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interest as this would indicate multiple solutions. Multi- 
ple solutions are excluded from the problems prsented here- 
in since they are derived from Maxwell's equations and 
appropriate boundary conditions and, as such, uniqueness is 
guaranteed. The case that the denominator of equation (68) 
equals to zero requires closer examination. This case 
would yield an a •*■«> . This implies that x R ^ 0 as seen from 
equation (63) and this case is a trivial one for the pro- 
blems of interest. Therefore, this special case is not con 
sidered. The a and 8 determined above then allow the 

iteration equation (63) to be solved for x . = F(x ,y ). 

n+1 n J n 

With this corrected value of x + . , a similar set of con- 
traction correctors may be calculated for y n and y n+1 

solved for. The contraction coefficients y and 8 are cal- 
culated as: 


. I * (1 ' Y> Sy (x n+l’V 

V’WV 


for f y <5 WV * 0 <71) 


and 


g y (x n+l»y n ) 

lf * ( ’ S ' +1 ’ y ” > ' 11 Vx n+ 1 .y n > ' B * <X ” +1,yn> 


Y = 


tg/x n+ i.y n )- 1 n fy ( Sn+ 1 ,y n )-l] 

VVt-l'V 


" B x <x n+l , yn > 


(72) 
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When f y (X n+l’ y n ) “ °’ arguments similar to those previous 
hold and the following correctors are obtained as: 


Y 


SyW+l-V 

8n (x n + l- y n ) * 1 


(73) 


3 


(y -1) 6 x <* n +i-y n > 

1 ' f x U n.l' y n 


(74) 


The y and 3 are then used in equation (64) to find y n+1 m 

G(x n+ ^,y n ). The iterative process is then repeated until a 
suitable accuracy has been achieved. An example of this 
two-dimensional contraction corrector scheme is laid out in 

Table 4. 


Even though the foregoing text and results of Table 4 

0 

have indicated that the contraction corrector method may be 
successfully applied to analytical two variable problems, 
it remains to be seen that the theory is directly extend- 
able to numerical methods. Current results of the applica- 
tion of equations (63) and (64) with the correctors of equa- 
tions (67) through (74) as applied to the one-dimensional 
problems of Sections IV and V agree with the previous data. 
However, the application of the contraction corrector 
scheme to a two-dimensional canonical problem has yielded 
poor results. The problem appears to be the numerical evalu- 
ation of the partial derivatives. Without the accurate de- 
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f(x,y)=x= 2 x+ 2 y g(x,y )=y=x 2 + 2 y 

F(x,y)=x=ax+(l-a)f (x,y)+ 6 y- 6 g(x,y) 
G(x,y)=y=yy+(l-y)g(x,y)+3x-3f(x,y) 


n 

x n 

y n 


6< X „'V 

0 

1 

1 

4 

3 

1 

4 

3 

14 

22 

2 

14 

22 

72 

240 


3 diverges 


n x n y n a g Y 3 G(x n> y n ) 

01 1 0.67 0.67 0.5 

1 0.5 1 - - 0 1 - -0.25 

2 0.5 -0.25 - 

* * 

x =0.5 y =-0.25 


Table 4. Example of a two-variable contraction 
corrector method. 




68 


termination of f , f , g , and g , it is doubtful the con- 

x y x y 

traction corrector method will be successful in arriving at 
a convergent solution. Another problem that is not readily 
apparent is the stringent conditions placed on the partial 
derivatives. Invoking such strong conditions as given in 
equation (60) may cause computational problems and eased or 
adjusted conditions may be necessary uo reach a viable so- 
lution method. These problems are being investigated and 
useful results are expected in the near term [25] . 

The usefulness of the two-dimensinal contraction cor- 
rector can not be underestimated. For instance, when a 
plane wave strikes a one-dimensional grating at an inci- 
dence angle that induces current along two axes, a cross- 
polarized scattered field will be generated. The two- 
dimensional nature of the problem must be addressed to ob- 
tain the valuable cross-polarized information. The solu- 
tion for reflection of waves from actual mesh surfaces is 
also of interest and involves the use of the method de- 
tailed in this dissertation. 

An illustration of an actual mesh surface is shown in 
Figure 20. The material is molybedenum wire coated with 
gold and woven in a double-knit fashion. Reflector sur- 
faces of space-deployable antennas have surfaces of this 
nature and their reflection characteristics must be known 
to accurately predict their electrical performance. Figure 
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Figure 20. 


Actual mesh surface. 
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21 illustrates the periodic nature of the mesh surface. It 
can be seen from the outlined boxes that Floquet cells may 
be formed. This periodic nature may be capitalized on and 
used to formulate the mesh surface into a suitable problem 
for the contraction corrector scheme. 

Figure 22 is the discrete representation of one single 
mesh path. The discretation of the mesh path is accomp- 
lished with twenty-eight, approximately equal length, 
straight segments. With these segments defined, it is pos- 
sible to overlay a grid onto the segments to achieve a rep- 
resentation of the mesh surface. The equally spaced points 
on the grid represent the sampled points required for the 

Fourier Transform of the solution method. This typical 
grid representation of the mesh surface will allow the de- 
termination of the reflection characteristics of the sur- 
face. Other surface weaves and materials may be placed on 
similar grids and their reflection characteristics calcu- 
lated. Figure 23 illustrates the application of the grid 
to the segmented single mesh path. 
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actual mesh surface. 




I 


Figure 23. Grid representation of the segmented 
single mesh path. 
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VII. CONCLUSIONS 

A new method called the contraction corrector method 
has been developed to insure the convergence of the one- 
dimensional spectral-iteration approach when solving elec- 
tromagnetic scattering problems. The method was presented 
beginning with basic examples of the iterative method and 
progressing to detailed iterative operator theory. Comput- 
ed data generated from the scheme was critiqued and compar- 
ed to work previously presented and was found to agree well. 
The ability to predict the regions of solvability with the 
contraction corrector scheme was demonstrated. The method 
was then utilized to solve for the reflection coefficients 
from an infinite grating of thin wires and the results were 
compared to previously published data. The ideas of gener- 
alization to a two-dimensional problem were discussed. A 
solution technique for the two variable contraction cor- 
rector method was developed and numerical difficulties in 
applying the method were discussed. A sketch of an actual 
mesh surface was presented and the application of the spec- 
tral-iteration approach was detailed. A general reference 
list and an over-view of both the k-space problem formula- 
tion and the spectral-iteration approach were also included. 

The material presented in this dissertation gives a 
solid foundation for future research, and indicates that a 
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suitable method for analyzing scattering from arbitrary 
mesh surfaces can be developed. The recommendations for 
future research should be directed toward understanding the 
the shortcomings of the two-dimensional problem as detailed 
in Chapter VI. With the proper constraints placed on the 
two-dimensional problem and the correct evaluation of the 
partial derivatives, a suitable method of solution is deemed 
obtainable. Another alternative is the reformulation of 
the problem with another method of attack. Many numerical 
methods have yet to be applied to the problem and one may 
will be suited for the two-dimensional scattering problem. 
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The FORTRAN source code that follows will allow the 
user to solve for the electric field in the aperture of 
an infinite grating of thin wires as well as the reflection 
coefficient from the grating. The program is written in 
FORTRAN 4 and should run on most FORTRAN compilers with no 
modifications except the unit numbers on the read and write 
statements. The plot routine must be user supplied and the 
program modified accordingly. The input variables are spe- 
cified in comments internal to the program. Comments are 
also included to help the user follow the flow of the pro- 
gram. Some of the variables are self explanitory and the 
others are listed below. A suitable DEL, the derivative 
increment, has been found to be (0.01,0.01) for most cases. 
E(I) Electric field in the aperture 

JC(I) Current density on the strip 

G(I) Transformed Green function 

EUT(I) Dummy E field variable 

CREF Complex reflection coefficient 

HI Incident H field 

ED(I) Incremented E field in the aperture 
ALF(I) Contraction corrector 

CK? Constants used in the iterative equation 

K Propagation constant 

EOUT(I) Output variable used for the external plot 

Z Internal impedance of the strip 
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Lines 


Program Flow 


170-400 

410-520 

530-610 

620-640 

650-720 

740-920 

930-990 

1020-1110 

1120-1130 

1140-1150 

1160-1200 

1210-1230 

1260-1290 

1300-1370 

1380-1440 

1450-1510 

1520-1580 

1590-1610 

1630-1670 

1680-1740 

1750-1780 

1790-1840 

1850 

1880-1980 


Input data 

Initialize parameters 

Calculate incident E and H fields 

Calculate thin wire impedance 

Set up aperture/strip sample space 

Calculate Floquet propagation constants and 
Green function transform 

Initial aperture E estimate 

Si 

Calculate constants in iterative equation 
Store untransformed E 

a 

Take FFT of E a and (E & + DEL) 

Perform G FI 

a 

Take inverse FFT of G E 

a 

Perform complement truncation T (F[G E ]) 

C a 

Perform FFT on T (F[G E]) 

C 2L 

Perform G" 1 (F[T (F[G E ])]) 

C ct 

Calculate reflection coefficient 

Perform inverse FFT on G -1 (F[T (F[G El)]) 

Plot set up 

Calculate contraction corrector 
Apply contraction corrector if desired 
Calculate finite conductivity effects 
Perform FFT on first iterated E 

a 

Test for five iterations, next iteration 
go to 1160 

Output external plot and stop execution. 



0REG3WAL PAGE S3 

OF POOR QUALITY 

c 

£«***«* FINAL. S ♦*** DERIVATIVE CORRECT ****** 
C Author Jerry C. Brznri 

C DATE August 12, 19 c '3 

C All richts reserved 
C 
C 

C Dimension all arrays 

COMPLEX E(51?),JCC512>,&(512) ,EUT(512) ,DEL 
COMPLEX RI,CPEF,HI 

COMPLEX ED(512),ALF(51?),FRIM£(512>,UK512>,EDUT(5 4 2> 
COMPLEX J,CKl (512) ,CK 2(512), CKC512) ,7 
DIMENSION POUT (512 ,10,2 ) 

REAL K,K2 

C a = r loquet cell dimension 
C F = Strip size 

WRITE (6,*) * INPUT FLOQLET CELL SIZE , STRIP SIZF' 

WRITE (6,*) ' NORMALIZED IN WAVELENGTHS' 

READ (5 ,* )A,B 

C FREQ a Frequency in Hertz 

WRITEC6,*) ' INPUT FREQUENCY IN HERTZ' 

RcAD(5,*) FP£<? 

C MAX = F FT size = Number of saircles per cell 
C Iw = loy2 (max ) ; i.e. FAX = 2**li. 

C IP = D no plct; IP = 1 p Ic t; 

C NOTE External plot routine reouired 

C I A C = C No correction; 

C IAP«1 Contraction corrector scheme applied 

WRITE (6,*) 'INPUT h, LOGZCO, PLOT OPTION, ADJUST OPTION' 
READC5,*) MAX, IU ,1 F ,1 AD 
C SIP * Conductance cf strip 

WRITEC6,*) ' INPUT CONDUCTANCE OF STRIP IN SIEMANS' 
READ15,*) S1G 

C TH = Theta angle of incidence 

C DEL = Complex increment usee ir derivative routine 
C PH = Phi angle o* Incicence 

WRITE (£,*) ' INPUT TH£TA ANGLE, DERIVATIVE DELTA, PHI ANGLE 
48? R£AD(5,*)TH , DEL, PH 
WRITE(6, 1 6) PEL 

16 FORMAT!'-',' DEL= ',iE10.4> 

IF(IAD.GE.I) WRITE (6,299) 

29« F0RMAT('-',' ADJUSTMENT USED ') 

C Initialize routine constants 
PI=3.141 59T 
TPI=PI*2. 

C*2.997956E+8 
UU=4.E-7*PI 
EP s 8.854E-z 2 
ET* s SQRT (UU/EP) 

I T E R = C 

J= ( u. G , 1 .0? 

AL«M5=C/ FREQ 
RD=1?C./°I 
DR=1 . /«D 

C Calculate incident electric and magnetic field components 
C for electric field parallel to vires, i.e. no cross- 
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C 


CO tar i 2 at ion includec 
I F CPH .LT .45 • ) EINC=1.C 
IFCPH .6E.45 .) EINC=CCSC TH»C«) 
I F CPH .LT • 4* • ) STH = C. 


ORIGINAL PAGE is 

Of POOR QUALITY 


I F CPH *6T .45 • ) STH = TH 
H1*1./ETA*CC0SCTH* CR)+J*SIMTH*DR) ) 

IFCPH.LT. 45.) h:=1./FTA*C0S(Th*DP) 
Z=SQPT(TPI*FRE0*UU/2./SlG)'*(1. r ' t 1.0)/TPI/b 
WRITE (£ .15) Z 

15 FORMAT C ' Z=C,F10.? .E10.2 

Calculate number of samples on strip and in aperture 
t All = A -P 

NsIFlXCTAU/A+FLOATCKf X)> 


N1=N+1 


NMl=N-1 


MAX1 = MX*1 


UR I T E ( 6 . 3 3) * . T . T A l! .FREL.TH.N , * A X 

TO F0RMAT('-',3F10.5,E1^.?,lFli.4 f 2nC) 
I FCN1 .GT.MAX) GOTO 9^8 


k=tpi/alakp 

K 2 =K*K 

SK=K*SINCTH*PR)*COS ( n H * L P ) 
SSf = X’*SI\C' r H*D<’)*SIN(PH«pR) 


W=TPI*FREQ 

maxi=max-1 

C Calculate Green function transform 
DO 4 C 1=1, M A X 
IFCI. GT. MAX/2+1 )GOTO 5? 

U=TPI*(I-1)/A-SK 
GOTO 63 

5 0 U = TPI«(I-MAX-1> /A-SX 

6 J U = U*U-tSSX*SSK 

IFCU.GE .X2) GOTO 7£ 

G (I )=-J*5QRT(K2-U) 

GOTO 44 

70 G( I) =-SQRTCU-K2 ) 

44 6<1)*G(X)-SSK*SSK/E(!> 

4 C CONTINUE 

C Initial F field estimate 
DO 32C 1=1 ,N 
E(I)=C1.C,r.G) 

320 EDCI)*E<I)*DFL 
CO 33 C I = N 1 i* AX 
EDC> = CG.Q,G.C) 

330 E (I )= (0.0 ,?. n ) 

C NOTE...* Iterative fore usee in this program is x = AX+E 
C Calculate M 8" ccrtion cf iterative equation 
DO 11C 1=1. MAX 
11C CK 1 (I ) = HI*J*VI*LU 

CALL FFTCCX1.IW.VAX) 

DO 120 I=N1,"AX 
1 2 C CK? Cl )=HI *W*UU/ J 
DO 14£ I =1 ,N 
1 4 C CK?(I)=<O.C,P.O) 

CALL FFTCCXi.IW.MAX) 

DO 1 3 C 1=1, MAX 

13C CKCI) = CCX1CD + CK2C1))/GII) 
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00 30C 1=1, MAX 
3DC UMI)=E(J) 

CALL f FT (£0 ,IW ,MA X J 
CALL F FT(E, IW ,MAX) 

C Begin "AX" portion of iterative equation 
C Perform (6*E) 

80 00 10C 1 = 1, MAX 

E0U) = C0*UG<PDCI>*C(I )> 

IOC EII)=CCNJG(E( I)*G( I)) 

C Perfor® inverse transfer* of (G«E) 

CALL FFT(E,IW,MAX) 

CALL FFT CEO ,IW ,MAX J 
WRITE (3,1*59) ITER 
199 FORMAT(1CX,110) 

C Perfor® truncation operation T(G*£) 

00 150 I =1 ,N 
ED ( I ) = (0 • 0, C • 0 ) 

15C E(I)=(C. 0,0.0) 

DO 16 C 1 =N 1 »K AX 
ED (I ) =C ON JG (ED (I)) /M^X 
E(I)=CCNJG(E( I))/M/X 
C Calculate current density on strip 
1 6 C JC(I) = EC!)*J/W/UU-H 
C Perform inverse transfer* cn T(6*t) 

CALL FFT(£, IW ,*AX) 

CALL FFT(ED,IW,MAX) 

C Perform T(G*E)/G and acd constant ,r O M 
DO 17 C 1 = 1, MAX 
E (! > = E (I )/G (I ) + CK ( ] ) 

ED<I>=ED(I)/G(I)*CK(I) 

EDUT( 1) = C0NJG<ED(1 )>/«AX 
EUT(I)=CCNJG(E<I)> / MAX 
17C CONTINUE 

C Calculate reflection coefficient 
°I=J*SIN(S T H*DR)/CCS(STH*LP) 

CREF=(E(1)/KAx*ElN()^J*S!L(STH*DR)*APS(1.-ABS<E(1)/MAX*tIf4C)) 
C«tF = CPEP/(CCS(STH<D: )* J * S IN ( $TH * D R ) ) 

REP=CAP$(C°EF) 

WR ITE (6 , A 5? ) RE F , C F E F 

453 FORMAT(6X,' REF= * , Fl u . i , £0X , 'CRE F= ' , FI 0. 2 ,2 X , F 1 C . i ) 

C Perform inverse t r ans f c r* at ion to obtain first iterated 
C electric field 

CALL FFTCEUT, IU,MA>) 

CALL FFT (EDU T »IW ,M AX) 

DO 2CC I =1 , M A X 
EDUT(I)=C0NJG(EDUT(1)) 

20C EUT(I) = CONlJG(EUT(I)) 

C Plct option 

XfClP.LT.05 GOTO 150 

CALL JPPLOT (EUT ,MA),TTEh,fcCuT) 

180 ITER=IT£P+1 

C Call subroutine to calculate derivative of operator and 
C contraction factor 

CALL DERIV(EUT,EDUt,"AX , D R IN E , 4L F , I A D , D EL , C ON ) 

WRITE (6,340) CON 

34 C FOPHATC/'IOX* * CONTRACTION = *,Flu.4) 

C Perfcr* correction on electric field; NOTE... If I A D = 0 
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ORIGINAL PAGE" IS 

no adjustment is usee as ALF(1)=3.0 ^ POOR QUALITY 

DO 35 C 1 = 1, N 

EUT<n*ALFfI>*UlCl)-»C1.-»LF(I))*EUTCI) 
tl1(I) = EUT(I) 

350 EDUT(I)=EU T (I ) + D£L 
310 CONTINUE 

C Calculate field on strip due to finite conductivity 
DO 21C I*N1,MAX 
EDIIT(I)=“JC(I)*Z*B 
21 C EUT(I) = -JCf I)*Z*B 

C Perform transform of E fielc tc begin next iteration 
CALL FFT(EUT,IW,MAX) 

CALL FFT CEDUT ,IW,MAX) 

DO 22G 1 = 1, FAX 
ED(I)=EDUTCI) 

22 C E(I) = El'T(I) 

IFC1TEP.6T.5I GOTO 24C 
GOTO £C 
240 COMTIME 

DO 777 1 = 1 , MAX 

PINDEX = FL0AT( 1-1 ) / FLOAT IMAX-1 )*A 
C Perform plot ftretiens fer external Plot routine 

WRITE (? 2, 77 6) PINDEX,EOiriI, 1,1) ,E0UT(I,2,1 ) ,£OUT(I,3, 1) 
4 t ECUT (1,4,1 ) ,£0UT ( 3 , f ,1 ) 

WRITE C23,77E1 RlNDEX,EOOT(x,1 f ?) , EOUTl 1,2,2) ,EOU T (1,3, 2) 
♦ ,ECUT<I f 4,7),ECUT( I , c , 2 ) 

77F F0PMATC6F8. 2) 

777 COMTIME 

WRITE 16,26 0) 

260 F0PKAT('-',1CX TIMLLY EXIT') 

59 E WRITE (6, 59) 

59 FORMATC'-',* ERROR IN N" ) 

5595 STCP 
END 

suproutime J«PL0T (A,!NT ,ITER,E0UT) 

COMPLEX A ( 1 ) 

DIMENSION AFP (512) , E OUT < 5 1 2 , 1 0 , 2 ) , PH A S E ( 5 1 2 ) 

N = ITE Ml 
DO 8C 1=1, INT 
AMP (1 ) = CABS (A (I ) ) 

PHY = A IFAG (A (I ) ) 

phx=real(a(i)) 

IF (AR S<PHX).LT.1 ,E-2: ) FEXs'i.t-Z" 1 
PHASE (I)=57«295*AT A 2 ( p H Y ,PHX ) 

EOUT ( I,N,1)=AKP(I) 

E 01! T ( I,N,2) =PH ASE ( I) 

80 CONTINUE 
RETURN 
END 

SUBROUTINE DERIV(E,ED, m AX,PkIHE,ALF,IAD,DEL,COM) 

C This subroutine calculates the derivative and 
C contraction factor ol the operator 

C0MPL EX E (MAX ) f ED (FAX) ,ALF (MA X) ,PRIMF (MAX) ,DEL 

C 0 N = C .0 

DO 1^ 1=1, MAX 

FRIME(I)=(ED(I)-E<l))/DfcL 

COM=CCN+CCAfiS (PRIMF (I)*DEL))*»? 


85 


/»LFCI) = P*?If^E< IJ/CPPI^EC 
13 XFCX«D*LT.1> ALFll)*C.O 

CON= SORT (CCS/ FLO AT (MAX) ) / C *bS (DEL) 

PETURN 

END 

SUBROU T IN£ FFTCAtMtN) 

C This Is the F FT subroutine called From the main orogram 
COMPLEX A(S) t U,W*T 
... N = 2 * * M 
NV2=N/2 
NMl = n -1 
J = 1 

00 7 1*1 * NM 1 
IFCI.6E.J) GOTO 5 
T*« (J) 

A(J)*ACI) 

A (I )=T 

5 K = NV2 

6 IF(K.GE.J) GOTO ? 

J = J - K 

K =K / 2 
GOTO 6 

7 J=J+K 

P1=3*1A159?6535&°7^3 

00 20 L*1.m 
LE=2**L 
LE 1 *L E*/2 
U*(1 ,C,u.O) 

W = CMPLX(COS (P 1/LEI ),0 IS (Pi/cFD) 

00 20 J * 1 » t. E 1 
DC 1C I*J.*ilE 
IP=I +L61 
Ts A (IP-) *U 
AIIP)*A< !>-T 
ID A(I)=A(I)+ T 
20 U=u*u 
RETURN 
END 
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APPENDIX B 


The general form of the basic iterative process is 
given as 



The linear operations indicated in equation (Bl) may be 
easier to understand if a flow of the operations is presented. 
Figure B-l illustrates the iterative process used in the 
computer program of Appendix A. The initial estimate of 
the electric field in the aperture is transformed and the 
transformed Green function is applied along with Floquet’s 
theorem. The inverse transform is applied to this result 
to obtain the current density on the strip. The inverse 
of the previous steps is then performed to obtain the 
first iteration of the electric field in the aperture. 

That is, the current density is truncated over the aperture 
and transformed, the inverse transform Green function is 
then applied along with Floquet's theorem to obtain the 
iterated aperture electric field. This iterated electric 
field may then be truncated and the process repeated until 
a suitable solution and tolerable error is .obtained. 


Figure B-2 illustrates where the basic iterative pro- 



H FOURIER TRANSFORM I I INVERSE II II 

AND APPLY I 1 FOURIER l^J OBTAIN ITERATED [^|»| FINAL 

FLOQUET'S THEOREM I I TRANSFORM I I CURRENT DENSITY II I RESULT 
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cess may be interupted to accomodate the contraction cor- 
rector scheme. This point is chosen since an iterated 
solution has been generated. With this iteration and any 
previous iterations a filtering process may be employed. 

The text pointed out the usefulness of the contraction 
corrector scheme since it guarantees a convergent solution. 
Thus, the iterative process can be viewed as a feedback 
process whereby the errors of processing can be eliminated 
and a useful solution of the problem obtained. 


C- - 
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APPENDIX C 

The optimum contraction corrector for the bne-dimen- 
sional problem is shown to be equivalent to Newton's 
method. The proof is as follows: 


f( V * 0 


G(x n> = x n+l * + (1 ' R) 9<V 

g’ g' 

= x n + 1_ 9 (x n> 

g 1 g'-l-g 1 

= Fi V + T'-i 9(x n ) 

g' ” g(x n ) 

g' -1 x n + g'-l 

, 9' x n - 

g'-l 

. 9’ x n - 9<* n > 4 x n - x n 

g'-l 

(g'-l) x x - g (x) 

— n .. n n 

(g'-l) g'-l 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6 ) 
(C7 ) 
(C8) 



